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Influence of quantum critical point (QCP) of first-order 
valence transition (FOVT) on Ce- and Yb-based heavy 
fermions is discussed as a key origin of anomalies such 
as non-Fermi liquid, metamagnetism, and unconven- 
tional superconductivity. Even in intermediate-valence 
materials, the QCP of the FOVT is shown to be induced 
by applying the magnetic field, which creates a new char- 
acteristic energy distinct from the Kondo temperature. 



It is stressed that the key concept is closeness to the QCP 
of the FOVT by pointing out that the proximity of the 
QCP explains sharp contrast between X=Ag and X=Cd 
in YbXCu4, field-induced valence crossover in X=Au, 
and field-induced first order transition as well as non- 
Fermi-liquid critical behaviours in Celrlns. The valence 
fluctuations can be a key origin of unresolved phenom- 
ena in this family of materials. 
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1 Introduction Quantum critical phenomena in itin- 
erant fermion systems have been intensively studied in the 
context of spin fluctuations [1-3]. Recently, quantum criti- 
cal phenomena in charge degrees of freedom have attracted 
much attention [4]. One of such highlights is valence tran- 
sition, which is a phase transition with a valence of ma- 
terials element showing a discontinuous jump. A typical 
example is known as the 7- a transition in Ce metal, where 
the isostructural first-order valence transition (FOVT) oc- 
curs in temperature-pressure (T, P) phase diagram [5] (see 
Fig. 1(a)). The valence of Ce changes discontinuously be- 
tween Ce+3 "3 phase) and Ce+^-'^'^ (a phase) at T = 300 
K [6]. The FOVT line terminates at the critical end point 
at (T, P) = (600 K, 2 GPa). As diverging density fluctu- 
ations at the critical end point in the liquid-gas transition, 
valence fluctuations diverge at the critical end point. In Ce 
metal, the critical-end temperature is so high that quantum 
criticality has not been discovered. 

When the critical-end temperature is suppressed to 
T = K by changing materials parameters, the quan- 
tum critical end point emerges. Since the quantum crit- 
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Figure 1 (color online) T-P phase diagram of (a) Ce and 
(b) CeCu2(Sii_a;Ge2;)2. (a) FOVT between the 7 and a 
phases terminates at the critical end point (CEP), (b) Uni- 
fied phase diagram of a: = [8], x = 0.1 [10], and 
X = 1 [7]. As for P axis, P^ = 4.5 GPa {x = 0), 
P^ = 5.5 GPa (x = 0.1), and P^ ^ 15.5 GPa {x = 1) 
are set as the same position (see ref. [10] for details). The 
dotted line at P^ is a guide for the eyes. 
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ical point (QCP) is defined as the point at which the 
continuous-transition temperature touches T = K, we 
refer to the quantum critical end point as the QCP in this 
paper The quantum criticality and electron instabilities 
arisen from the QCP of the FOVT has been recognized 
as emergence of two superconducting domes in the T- 
P phase diagram of CeCuaGca [7], CeCu2Si2 [8,9], and 
CeCu2(Si:EGei_^)2 [10,11]: Near the QCP P ^ Pc where 
the Neel temperature is suppressed under pressure, a su- 
perconducting phase emerges (SCI). Interestingly, further 
enhancement of the superconducting-transition tempera- 
ture Tsc emerges by further applying pressure (SC2) (see 
Fig. 1(b)). Around the pressure P — P^ where Tsc has 
a maximum in the SC2, the quantum critical behaviours 
such as a T-linear resistivity and a remarkable enhance- 
ment of residual resistivity have been observed. Since the 
resistivity shows the dependence near P = Pc, i.e., 
the criticality predicted theoretically for the antiferromag- 
netic (AF) QCP in three spatial dimension [1-3], these 
anomalies near P = P^ should be ascribed to a distinct 
origin of AF quantum criticality. 

Theoretically, a possibility of emergence of valence 
criticality at P = Pv has been pointed out in ref. [12]: 
A key experimental feature is the pressure dependence of 
the coefficient A estimated from the low-temperature 
part of resistivity [4]. A remarkable fact is that A sharply 
decreases with 2-3 orders of magnitude, when P in- 
creases and across P = P^ [7-11]. Since A is propor- 
tional to 7p with 7c being Sommerfeld constant by so- 
called Kadowaki-Woods relation [13], this suggests that 
the effective mass of electrons to* /mo is drastically re- 
duced at P = Pv This sharp crossover is also reflected 
in the sharp change of the Kadowaki-Woods ratio between 
strong-correlation regime A/j'^ = 10^'^ (P < Pv) and 
weak-correlation regime A/-j^ = 10^^ (P > Py). If 
we remind the fact that the effective mass is expressed 
by the /-electron number per site in Ce (and also Yb 
for the /-hole number) based heavy fermion systems: 
m* /nio oc (1 — n//2)/(l — nj) shown by the Gutzwiller 
arguments for periodic Anderson model [14], the above 
observations imply that the /-electron number is sharply 
deviated from Uf ^ 1 when P exceeds P^. In other words, 
the valence of Ce increases sharply from nearly trivalent 
3-1- (n/ ~ 1) realized in the low-pressure regime to a larger 
valence 3 + 6{nf^l — S) with S being a positive value. 

Recent ^■^Cu NQR measurement in CeCu2Si2 under 
pressure has detected that near P = Pv ~ 4.5 GPa the 
NQR frequency i^q starts to deviate from linear increase 
expected from the monotonic volume shrinkage under hy- 
drostatic pressure [9]. The downward deviation from the 
monotonic increase in the i'q-P plane indicates that the 
charge distribution around the Cu site changes around P = 
Pv, strongly suggesting that the Ce valence changes at 
P = Pv. Enhanced valence fluctuations have been shown 
theoretically to make the impurity scattering strong, giving 
rise to the enhancement of the residual resitivity [15]. The 



local nature of valence fluctuations has been also shown to 
cause the non-Fermi-liquid critical behaviour in the resis- 
tivity: The T-linear dependence emerges near the QCP of 
the FOVT [8]. 

In Fig. 1(b), Tsc in the SC2 starts to increase at the 
lower pressure than Pv, i.e., just before the Ce valence in- 
creases. After the sharp valence crossover for P > Pv, the 
superconductivity is immediately suppressed. This feature 
is reproduced by the slave-boson mean-field theory taking 
into account of Gaussian fluctuations applied to the peri- 
odic Anderson model with an inter-orbital Coulomb repul- 
sion (see Eq. (1) below) [16]: In the lower pressure region 
than Pv, Tsc is enhanced. The density matrix renormaliza- 
tion group (DMRG) calculation in the same model in one 
spatial dimension also shows that the enhancement of the 
superconducting correlation occurs for P < Pv near the 
QCP [17]. It is noted that the valence-fluctuation mediated 
pairing mechanism proposed in ref. [12] obtained the firm 
grounds by unbiased numerical calculation [17]. 

Yb systems provide a hole analog of Ce systems, where 
in the 4/ shell Yb'^^ gives 13 electrons and Yb^+ gives 
fully occupied 14 electrons. Thus, in the hole picture nj ~ 
1 {uf = 0) is reaUzed for Yh^+ {Yh^+). In Yb sys- 
tems, the isostructural FOVT has been also known to oc- 
cur in YbInCu4 [18,19]: As T decreases, the first-order 
transition from Yb+^-S^ (^nf ^ 0.97) to Yb+2-84 (^^^ ^ 
0.84) takes place at T = 42 K [20]. In the series of 
YbXCu4, when X=In is replaced to the other elements such 
as X=Ag [19,21] and Au [22], the first-order transition 
has not been observed. However, anomalous behaviours 
such as metamagnetism [19] and emergence of new char- 
acteristic energy scale distinct from the Kondo tempera- 
ture [19,22] have been observed, which seem to be related 
to the enhanced valence fluctuations. 

The critical behaviours and electronic instabilities ob- 
served in the Ce- and Yb-based systems seem to indicate 
a strong influence of the QCP of the valence transition. As 
seen in Fig. 1(a), the FOVT terminates at the critical end 
point in the T-P phase diagram in Ce metal. When the 
critical-end temperature is suppressed by controlling ma- 
terials parameters, and enters into the Fermi degeneracy 
regime, the diverging valence fluctuations are considered 
to be coupled to the Fermi-surface instability. This com- 
bined effect of critical fluctuations and quantum fluctua- 
tions seems to be a key mechanism in understanding the 
anomalous Ce- and Yb-based heavy fermions. In the fol- 
lowing sections, we show how the QCP of valence tran- 
sition is controlled by the magnetic field, and discuss its 
influence on Ce- and Yb-based systems in the subsequent 
section. 

2 Field-induced QCP of valence transition We 

analyze the magnetic field dependence of the QCP of the 
FOVT on the basis of the extended periodic Anderson 
model as a simplest minimal model for Ce- and Yb-based 
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heavy fermions: 



1.50 



H ^ Hr + Hi 



hyb 



Ufa 



sn (1) 



where He = 

The Ufc term is the Coulomb repulsion between / and 
conduction electrons, and is considered to play an impor- 
tant role in the valence transition. For example, in the case 
of Ce metal which exhibits the j-a transition, the 4/- and 
5(i-electron bands are located at the Fermi level [23]. Since 
both the orbitals are located on the same Ce site, this term 
cannot be neglected. In the case of YblnCu4, Hjjf^ also 
plays a crucial role for the FOVT in the hole picture of 
Eq. (1) [24]. Most of Ce- and Yb-based compounds seem 
to have moderate values of U fc- However, they seem to be 
affected by close proximity to the QCP of the FOVT as 
discussed below, which requires a moderate magnitude of 
Ufc in Eq. (1). 

By applying the slave-boson mean-field theory [16] 
to Eq. (1), we have determined the locus of QCP of the 
FOVT under the magnetic field [25]. Figure 2 shows the 
ground-state phase diagram in the e f-Ufc plane for D — 1, 
V ~ 0.5, and C/ = oo at n = {uf + nc)/2 = 7/8, where 
the conduction band is set as = fc^/(2m) — I? in three 
spatial dimension, and ric is a conduction-electron number 
per site. At h = 0, when is deep enough, the Kondo state 
with Uf = lis realized. When e / approaches the Fermi 
level, / electrons are moved into the conduction band via 
hybridization, giving rise to the mixed-valence state with 
Uf < 1. The FOVT between them is caused by Ufc, since 
large Ufc forces electrons to pour from the / level into the 
conduction band. The FOVT line (solid line with open tri- 
angles) terminates at the QCP (filled circle), where the va- 
lence susceptibility Xv = ~dnf / de f, i.e., valence fluctua- 
tion, diverges. The dashed line represents the points where 
Xv has a maximum as a function of ef for each U fc, im- 
plying that even in the valence-crossover regime valence 
fluctuations are well developed. Note that the hybridiza- 
tion between f and conduction electrons is always finite in 
the Ef-Ufc plane shown in Fig. 2. Hence, the Fermi surface 
is always large [|T7]|. 

Under the magnetic field, the QCP shows non- 
monotonic field dependence as shown in Fig. 2: The QCP 
shifts to the smaller-[/ fc and larger-e/ direction for small 
h, and upturn behaviour emerges for h > 0.04. Since 
the Kondo temperature at the QCP (sj*^^ 



1.40 



QCP 



0.074, 
upturn 



(0.356, 1.464) under h^Ois estimated as Tj 
the characteristic field of the beginning of the 
roughly corresponds to T^'~^^. The origin of the upturn 
is due to the broadening of the spectral weight of /- and 
conduction-electrons Green function by the magnetic field 
larger than T^'~^^ [25]. If Ce- and Yb-based materials are 
located just on the locus of QCP's for h > (a gray line) 
in Fig. 2, the metamagnetism will be observed at h = 



1.30 



1.20 




Figure 2 (color online) Ground-state phase diagram in the 
e/-t//c plane for £» = l,V = 0.5, and [/ = ooatn 7/8. 
The FOVT Une with a QCP for h = 0.00 (open triangle), 
h = 0.01 (filled triangle) h — 0.02 (filled inverse triangle), 
h = 0.03 (filled star), h = 0.04 (filled diamond), h = 0.05 
(filled square) and h = 0.06 (open square). The gray line 
connects the QCP's under h, which is a guide for the eyes. 
The dashed lines represent the valence-crossover points at 
which Xv has a maximum as a function of for each Ufc- 
Inset: m-h curve for {sf.Ufc) = (-0.354,1.458) (thin 
line) and (-0.349,1.442) (bold line), m is defined by m = 



(as seen in inset of Fig. 2). It is remarked that at the field- 
induced valence QCP, not only valence susceptibility Xv, 
but also the magnetic susceptibility Xs = —dm/dh di- 
verges. Thus, uniform spin fluctuations at fc = are also 
enhanced. 

It is noted that the field-induced QCP of the FOVT has 
also been confirmed by the DMRG calculation in Eq. ( 1) in 
one dimension [25]. This indicates that not only unbiased 
calculation reproduces the above result, but also a differ- 
ence of spatial dimensions does not alter the main conclu- 
sion. This is ascribed to the locality of the valence transi- 
tion, which has a local atomic origin [26]. 

The metamagnetic field h = h^^ corresponds to the 
difference of Tk at the QCP between /i = and h = h^^: 



T^'^^ih ^ 0) - T^'^''{h = 0). This indicates 



hQCP 



that there appears a distinct energy scale from the Kondo 
temperature, which is characterized by the closeness to the 
QCP of the FOVT. As mentioned above, Ce metal is con- 
sidered to have a considerable value of Ufc due to its onsite 

origin. Hence, it seems to be located at Ufc > U'f'^^ in 
Fig. 2, giving rise to the FOVT. On the other hand, in Ce- 
based compounds, the energy bands located at the Fermi 
level usually consist of 4/ orbitals at the Ce site and non- 
4/ orbitals at the non-Ce site. For example, 4/ orbitals at 
the Ce site and 3p (4p) orbitals at the X site are respon- 
sible for the bands in CeCuaXz (X=Si (Ge)). Then, Ufc 
has the intersite origin in Ce compounds. Hence, most of 
Ce- and also Yb-based compounds are considered to be 
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located in the intermediate valence-crossover regime for 
Ufc < U'j^^^ in Fig. 2. This is consistent with the exper- 
imental fact that most of those compounds do not show 
the FOVT, but merely show the valence crossover How- 
ever, the results revealed in Fig. 2 indicate that even such 
compounds can be affected by the proximity to the valence 
QCP by applying the magnetic field. In the next section, we 
discuss how this newly clarified mechanism resolves the 
outstanding puzzles measured in Ce- and Yb-based com- 
pounds so far. 

3 Comparison with experiments 

3.1 YbAgCu4 and YbCdCu4 YbAgCu4 and YbCdCu4 
have paramagnetic-metal ground states. The analysis of 
T dependence of uniform magnetic susceptibilities of 
YbAgCu4 and YbCdCu4 concluded that both have nearly 
the same Kondo temperatures, about 200 K [19]. However, 
only in YbAgCu4 a peak structure appears in the mag- 
netic susceptibility around T ^ 40 K, but in YbCdCu4 
a conventional Pauli paramagnetism appears [19]. Such a 
contrast has been also detected in the magnetic-field re- 
sponse: YbAgCu4 shows a metamagnetic increase of the 
magnetization around /i ^ 40 T, while YbCdCu4 does 
not [19]. These observations suggest that there exist a 
distinct energy scale from the Kondo temperature. 

Our results explain this sharp contrast. Figure 3(a) 
shows the schematic contour plot of the /-hole number per 
site, Uf, which can be also regarded as the contour plot of 
the Kondo temperature, Tk oc (1 — — ?^f/2) [14]. 

In the small (large) e / and Ufc regime, the Kondo (mixed 
valence) state with Uf ~ 1 (nj < 1), i.e., small (large) 
Tk is realized. Since YbAgCu4 and YbCdCu4 have nearly 
the same Tk, both are considered to be located in the same 
contour area (see Fig. 3). If YbAgCu4 is closer to the QCP 
with a distance about /i ~ 40 T than YbCdCu4 which 
is less close to the QCP, the metamagnetic behaviour in 
YbAgCu4 at ft, 40 T can be naturally understood. 

This picture also gives a natural explanation for the ori- 
gin of the peak structure observed in the uniform suscepti- 
bility at T - 40 K in YbAgCu4. Figure 3(b) illustrates a 
schematic phase diagram of the T-e f-Ufc space. As shown 
in inset of Fig. 2, the metamagnetism with diverging mag- 
netic susceptibility emerges at the QCP as well as the crit- 
ical end points represented by the thick solid line with a 
circle in Fig. 3(b). An important result is that even at the 
valence-crossover surface extended from the FOVT sur- 
face in Fig. 3(b), the magnetic susceptibility Xs as well as 
the valence susceptibility Xv is enhanced [26]. If YbAgCu4 
touches the valence-crossover surface at about T* ^ 40 K, 
the observed peak in Xs{T) at T ~ 40 K is naturally ex- 
plained. The reason of no conspicuous peak in XsiT) in 
YbCdCu4 can be understood as too large temperature inter- 
val to touch the valence-crossover surface (see Fig. 3(b)). 
Indeed, a volume expansion was observed below T ~ 40 K 
in YbAgCu4 [21], which is directly related to the sharp 




Figure 3 (color online) (a) Schematic contour plot of rif, 
i.e., the Kondo temperature Tk in the e f-U fc plane for the 
model (!) at h ~ 0. The FOVT line (solid line) terminates 
at the QCP (filled circle). The QCP reaches YbAgCu4 by 
applying h about 40 T represented by an arrow. YbCdCu4 
is located far away from the QCP. The valence crossover 
surface reaches YbAuCu4 by applying h about 1.3 T (see 
text), (b) Schematic T-e f -Ufc phase diagram at ft, = 0. 
YbAgCu4 touches the valence-crossover surface at about 
T = 40 K. YbCdCu4 is located too far from the QCP 
so that it needs too large T interval to reach the valence- 
crossover surface. 



valence crossover from Yb+^-^f (T > 40 K) to Yb+^-S'^ 
(T < 40 K) [27]. 

3.2 YbAuCU4 Recent ''^Cu NQR measurement in 
YbAuCu4 has revealed that the characteristic temperature 
T* emerges in the T-h phase diagram as illustrated in 
Fig. 4(a) [22]. Under the magnetic field for ft > 1 T, as T 
decreases, ^'^Cu NQR frequency i^q decreases sharply at 
T*{h). This indicates that the electric-field gradient at the 
Cu site changes sharply at T*{h), implying that the charge 
distribution changes drastically around the Cu site. This 
strongly suggests that the valence of Yb changes at T*(ft). 
As ft increases and across ft* ^ 1.3 T, the coefficient 
A of the resistivity sharply decreases, and the residual re- 
sistivity has a cusp-like peak structure at ft* [22]. These 
observations also suggest that the Yb-valence crossover 
temperature T* is induced by applying the magnetic field. 

This can be naturally understood if YbAuCu4 is lo- 
cated just below the locus of the QCP as illustrated in 
Fig. 3(a). As shown in Fig. 2, for ft < T^*^^(ft ^ 0), 
the valence-crossover line extends to the Kondo regime 
(dashed lines for ft ~ 0.02 and 0.03 cross that for ft = 0). 
Then, by applying ft, the valence-crossover surface ap- 
proaches and reaches the point indicated as Au in Fig. 3(a) 
at ft = 1.3 T, which gives ft* in Fig. 4(a). As ft further 
increases, the QCP moves as an arrow in Fig. 3(a). Hence, 
the valence-crossover temperature T*(ft) increases as il- 
lustrated in Fig. 4(a). In the large ft limit, the Yb'^^ state 
is expected to be realized, since the Zeeman energy gain 
is mostly earned by the rif ~ 1 state. Hence, T*{h) is ex- 
pected to start to decrease for larger ft, and finally to touch 
absolute zero. Indeed, this behaviour should appear if the 
valence QCP's exhibit the upturn as shown in Fig. 2. 
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Figure 4 (color online) (a) T-Zi phase diagram of 
YbAuCu4 [22]. The Neel temperature Tn is suppressed 
and the Yb valence crossover temperature T* is induced 
for h > h* ^ 1.3T. (b) T-h phase diagram of Celrlng 
{h II c) [33]. The first-order transition is induced for h > 
hv (solid line). In (a) and (b). For T < Tfl (dashed line) 
the resistivity shows the Fermi liquid behaviour ~ T^. In 
(b) dotted lines are guides for the eyes. 



Here, we point out the close similarity of the T-h phase 
diagram between YbAuCu4 [22] and YbRh2Si2 [28]. In 
YbRh2Si2, a very similar phase diagram as in Fig. 4(a) 
has been obtained with remarkable anomalies near h ^ 
such as the T-linear resistivity, residual resistivity peak, 
changes of the Hall coefficient, magnetostriction, and 
magnetization [28]. We stress here that these properties 
are explained if the system is close to the valence QCP, 
as mentioned above. Indeed, a band-structure calculation 
has shown that the change of the Hall coefficient around 
h ^ hy m YbRh2Si2 can be explained by a tiny Yb 
valence change [29]. Hence, it is important to examine 
the possibility experimentally whether the Yb valence 
changes at in YbRh2Si2. Enhanced valence fluctuation 
is a possible origin of unconventional criticality observed 
in YbRh2Si2 [26]. 

3.3 Celrln5 Celrlns is a heavy fermion metal, which 
exhibits the superconducting transition at T ^ 0.4 K. 
In the T-h phase diagram, the first-order-transition Une 
emerges as illustrated in Fig. 4(b), which was detected by 
the jump in the magnetization curve [30-33]. Capan et al. 
observed anomalous behaviours that the residual resistiv- 
ity increases as h approaches hy. Furthermore, the power 
of the T dependence of resistivity shows a convex curve in 
the T^-^ plot, which seems to indicate that the resistivity 
tends to show T-linear dependence. They pointed out that 
non-Fermi liquid behaviours even evident in = may be 
related to the matamagnetic anomaly at ^ h^, although 
its mechanism was not clarified [33]. 

Our results give a natural explanation for these anoma- 
lies: This is readily understood if Celrln,5 is located in- 
side the enclosed area of the QCP line for 7^ as in 
Fig. 2. Namely, ?Xh — the system is considered to be lo- 
cated at the valence-crossover regime (i.e., for U < uf^^ 
in Fig. 2, since no evidence of the first-order transition 
was observed in any physical quantities as a function of 
T at h = 0.) However, when h is applied, the QCP of 



the FOVT approaches and eventually goes across, caus- 
ing the metamagnetic transition in the magnetization curve. 
As mentioned above, near the QCP the residual resistivity 
is enhanced [15] and the T-linear resistivity is observed 
in the wide-T region [8], which are quite consistent with 
the observations [30-33]. Furthermore, the emergence of 
the first-order transition line in the T-h phase diagram is 
also in agreement with the present picture. We stress that 
CeCoIns which has almost the same Fermi surfaces and 
similar crystalline-electric-field levels as those of Celrlns 
does not show such anomalies. This indicates that the view- 
point of the closeness to the valence QCP is also indispens- 
able in understanding the Cel 15 systems. 

4 Summary We clarified the mechanism of how the 
QCP of the FOVT is controlled by the magnetic field. 
Even in the intermediate valence materials, the QCP of 
the FOVT is shown to be induced by applying the mag- 
netic field. Valence fluctuations developed at the valence- 
crossover surface near the QCP cause various anoma- 
lies. We discussed that the proximity of the valence QCP 
resolves the outstanding puzzles observed in YbXCu4 
(X=Ag, Cd, Au) and Celrlns. A key concept is the close- 
ness to the QCP of the FOVT. Influence of the valence QCP 
offers a key origin in understanding unresolved anomalies 
in Ce- and Yb-based heavy fermions. 
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